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Generic First Order Orientation Transition of Vortex Lattices in Type II 

Superconductors 
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First order transition of vortex lattices (VL) observed in various superconductors with four- 
fold symmetry is explained microscopically by quasi-classical Eilenberger theory combined with 
non-local London theory. This transition is intrinsic in the generic successive VL phase transi- 
tion due to either gap or Fermi velocity anisotropics . This is also suggested by the electronic 
states around vortices. Ultimate origin of this phenomenon is attributed to some what hid- 
den frustrations of a spontaneous symmetry broken hexagonal VL on the underlying four-fold 
crystalline symmetry. 
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Morphology of vortex lattices (VL) in the Shubnikov 
(mixed) state of type II superconductors is not com- 
pletely understood microscopically. 1,2 ) Thorough under- 
standing of the VL symmetries and its orientation with 
respect to the underlying crystalline lattice are important 
both from the fundamental physics point of view because 
through those studies we can obtain the information on 
pairing symmetry (see below). It is also important from 
technological application of a superconductor, such as 
vortex pinning mechanism which ultimately determines 
critical current density. 

The first neutron diffraction experiment 3 -* was con- 
ducted for Nb to observe a periodic array of the 
Abrikosov vortices. Concurrently since then various ex- 
perimental methods have been developed, such as Bit- 
ter decoration of magnetic fine particles, /zSR, NMR or 
scanning tunneling microscope (STM). 1,2 ) The former 
ones are to observe spatial distribution of the magnetic 
field in VL and the last one is to observe the electronic 
structure of VL. The variety of experiments compile a 
large amount of information on VL morphology. Thus 
we need to synthesize many fragmented pieces of infor- 
mation which has been not done yet. Here we are going 
to investigate a mystery known for a quite some time and 
to demonstrate that the present methodology is powerful 
enough which might be useful for systematic studies of 
Shubnikov state in type II superconductors. 

As shown in Fig. 1 where the open angle 9 of a unit 
cell is plotted as a function of magnetic field B applied 
to the four-fold symmetric axis (001) in either cubic 
crystal (V 3 Si 4 < 5 )) or tetragonal crystals (YNi 2 B 2 C, 6 ~ 8 ) 
CeCoIn 5 9 '), the observed VL symmetries change in 
common, yet they belong to completely different su- 
perconducting classes (conventional or unconventional): 
Namely starting with the regular triangle lattice at the 
lower critical field B c \ , the VL transforms into a general 
hexagonal lattice with 60°, keeping the orientation 
a the same: parallel to (110). This direction is defined as 
(110) here, so that a || (110). The inset of Fig. 1 shows 
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Fig. 1. (Color online) Field dependence of open angle 6 for V3S1 
(Q) 5 4,5) YNi 2 B 2 C(CT) 6 - 8 ) and CeCoIn 5 (A). 9 > The open and 
filled symbols, respectively, belong to orientation a || (110) and 
a || (100). The shaded band indicates first order transition re- 
gion for these materials. Lines are theoretical curves by non-local 
London theory with /3 = 0.4 for guide of eyes. Inset shows a unit 
cell of VL with definitions of open angle 9, orientation direction 
a, and bonds A and B of VL. 



the definition of 9 and the orientation a. At a certain 
field B lst (shaded area in Fig. 1) whose values depends 
in external conditions, such as cooling rate, etc., the ori- 
entation a suddenly changes by angle 45° to a || (100) 
at a first order phase transition. Upon further increasing 
B, 9 continuously increases and finally becomes 9 = 90° 
at Bqq. This lock-in transition at B — Bgo is of sec- 
ond order. Those successive symmetry changes in VL are 
common for those compounds. Other A-15 compounds 
Nb 3 Sn 10 ) or TmNi 2 B 2 C 1:L) exhibit a similar trend. 

We explain this phenomenon by synthesizing two the- 
oretical frameworks: phenomenological non-local London 
theory 12-14 ) and microscopic Eilerberger theory. 15 ~ 17 ) 
The purpose of this study is to examine the universal 
behaviors of first order transition in the VL morphology 
under changing the degree of anisotropies coming from 
two sources: the Fermi velocity anisotropy and super- 
conducting gap anisotropy. In this analysis, we evaluate 
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the validity of non-local London theory, comparing with 
quantitative estimates by Eilenberger theory. Further, we 
discuss intrinsic reasons why the first order transition of 
VL orientation occurs, in viewpoints from frustrations 
and electronic states around vortices. Microscopic calcu- 
lation is necessary to discuss the electronic states. 

As we will see below, frustrations are a crucial key con- 
cept to understand this phenomenon. Antiferromagnetic 
spins on a triangular lattice is a well known example of 
frustration. Here a triangular VL is also a driving force 
due to hidden frustration which is more subtle than the 
spin case, leading to successive phase transition. Hexag- 
onal symmetry is incompatible with four-fold symmetry 
of underlying crystal lattice, because all bonds between 
vortices are not satisfied simultaneously to lower bonding 
energies. 

Quasi-classical Eilenberger theory 15-17 ) is valid for 
kp£ » 1 (hp: the Fermi wave number and £: the co- 
herence length) a condition met in the superconductors 
of interest here. Eilenberger equations read as 

K + «(*)■ (II + *A))/ = A(r,%, (1) 

K-t»(0)-(n-*A))/t = A*(r,% (2) 

with g 2 + fp = 1 for quasi-classical Green's func- 
tions g(u n ,r,8), f(u n ,r,6), and p(uj n ,r,9). In our for- 
mulation, we use Eilenberger unit where length, mag- 
netic field and temperature are, respectively, scaled by 
i?o = hvF /2nkBT C7 B — <p /2irR 2 , and transition tem- 
perature T c with flux quantum </> . Matsubara frequency 
Lo n = T(2n + 1) with integer n. The pairing interaction 
is assumed separable V(9,9') — Vo(p(6)cp(6') so that gap 
function is A(r,6) = ^{r)(f>{9). Since we consider two- 
dimensional case with cylindrical Fermi surface, we set 
the normalized Fermi velocity v(9) = vf/vfo as v {6) = 
v(#)(cos#, sin#). When we discuss the Fermi velocity 
anisotropy, we model it as v{9) = (1 + /3 cos A9) / y/\ — f3 2 
and (j>{6) = 1, which is called (3 model. 17 ) When we dis- 
cuss the gap anisotropy, v(9) = 1 and (f>{9) 2 = 1 — a cos 4$ 
as a-model (anisotropic s-wave), or (j)(9) = \/2cos2# as 
cZ-wave pairing. Here 9 is the polar angle relative to (100) 
axis. 

The selfconsistent equations for the gap function \I/(r) 
and vector-potential A(r) are 



*(r) = V N 2T ]T {4>*{0)f) 



w„>0 



V x V x A(r) 



~ J2 Im («s>- 



(3) 



(4) 



with V Q N = InT + 2T]T^ >0 w" 1 and Ginzburg Lan- 
dau parameter k. For average over Fermi surface, (...) = 
J Q 7r (. . .)d9/2nv(9) with extra factor l/v(9) coming from 
angle-resolved density of states N(9) = Nq/v(6) on 
Fermi surface. The self-consistent solution yields a com- 
plete set of the physical quantities: the spatial profiles of 
the order parameter ^(r) and the magnetic field B(r). 
The local density of sates (LDOS) for electronic states 
is calculated by N(r, E) = N a Re{g(Lj n ,r,9)\ iuJn ^ E+i0 ). 



The free energy density is given by 
F 



uj„>0 



K 2 B 2 (r)-Tj2 (^(A*/ + A/t) 



9, 



1 + 5 



Here, (. . .) is a spatial average within a unit cell of VL. 
Free energy should be minimized with respect to the VL 
symmetry and its orientation relative to the crystallo- 
graphic axes. In previous studies by Eilenberger theory, 
the first order transition of the VL orientation was not 
evaluated at low T, while only the transition from trian- 
gular to square VL was discussed. 16, 17 ^ In principle, we 
can investigate the whole space spanned by (B,T). But 
in practice it is not easy to exhaust the parameters in 
order to seek the desired physics. Thus our calculations 
are backed up by the non-local London theory. 

The non-local London theory is powerful and handy for 
B c i < B at low T. The nonlocal relation between current 
j and vector potential A in Fourier space is derived as 



L .ii{q) =i[qx A], 



: J2 m Qim(q)A m (q) (l,m = 



x,y) from the Eilenberger theory. 13 ) Since we use the 
Eilenberger unit here, the penetration depth A is changed 
to k in the length unit Rq (~ £). The kernel is 

\A(9)\ 2 v lVm 



Ql m (q) = 2TJ2 



Ul n >0 



(5) 



with ft 



A 2 (6), and A(6) is a uniform solution 



of the gap function. In previous studies by the non-local 
London theory for the first order transition in borocar- 
bides, higher order terms than q 4 was neglected. 12 ^ Here, 
we do not expand Qi m by q so that we can include all 
order contributions of q in Qi m - 1A ^ The corresponding 
London free energy is given by 

i 

+ {Q^)xxq 2 y - 2(Q- 1 ) xy q x q y }}-\ (6) 

where Q^ 1 is the inverse matrix of Qi m depending on q. 
The non-local London theory is valid near B c \ at lower 
T where the vortex core contribution is approximated as 
a cutoff parameter £ mimicking the finite core size effect. 
To explore wide ranges of (B,T) on a firm basis, we 
need to carefully check the validity of non-local London 
theory, using Eilenberger theory. We report our results 
for k = 89 and T = 0.2T C in both theories. The results 
do not depend on the k value unless k is approaching 
1/V2. 

First, we study the /3-model, i.e., anisotropy of Fermi 
velocity. Figure 2 shows a stereographic view of the open 
angle 9 as functions of B and (i. It is seen that (1) ir- 
respective of the (3 values the first order transition ex- 
ists seen as a jump of 9 in Fig. 2. (2) The jump of the 
first order transition becomes large as (3 increases. (3) 
(3 > /3 cr = 0.38 {(3 < f3 cr ) at T/T c = 0.2, the square lat- 
tice is (never) realized. (4) This square lattice becomes 
ultimately unstable for higher fields. Items (3) and (4) 
were also confirmed by Eilenberger calculation [see Fig. 
2(a) in Rcf. 17]. Note in passing that in CeCoIns the 
square lattice changes into a hexagonal lattice at a higher 
field, which reminds us the similarity, but we believe that 
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Fig. 2. (Color online) Stereographic view of the open angle 9 as 
functions of B and /3 in the /3-model by non-local London theory. 
First order transitions from 9 < 60° to 9 > 60° with 45° rotation 
of VL orientation are seen for any values of /3. 
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Fig. 4. (Color online) Zero energy LDOS N(r,E = 0)/N for 
hexagonal VL with 9 ~ 60° at B/Bgo = 0.07 in d x 2 _ H 2-wave 
pairing, where the gap minimum is along (110) direction, (a) 
Stable orientation a || (110). (b) Unstable orientation a || (100). 
Horizontal axis is parallel to a-direction. 




Fig. 3. (Color online) Open angle 9 as a function of B/Bgo for 
two orientations in the /3-model (/3 = 0.4, 0.5, 0.9) by non-Local 
London, and in the a-model (a = 0.5) and d-wave model by 
Eilenberger theory. Circles show 9 obtained by Eilenberger the- 
ory in the /3-model with /3 = 0.5. Arrows indicate £?i s t- Inset 
shows the closed packed square tiles and the associated hexago- 
nal unit cell, indicating 9 = 53° (lower inset) and 9 = 90° (upper 
inset) when B — » BgQ. 



it is caused by other effect, such as the Pauli paramag- 
netic effect. 19 ' 20 ^ 

When plotted as a function of B/Bgo by respective 
lock-in field Bgo as shown in Fig. 3, the ^-behaviors are 
similar for any /3 (> /3 cr ), while Bi st changes depending 
on (3. These universal behaviors well reproduce the ex- 
perimental data as shown in Fig. f. We also plot some 
points of 8 obtained by Eilenberger theory for (3 = 0.5 
(see the circle symbols in Fig. 3), which shows similar 
behavior as in the non-local London theory. Therefore, 
the non-local London theory can be applied reliably to 
superconductors whose anisotropy comes from the Fermi 
velocity. 

Next, we study the case of anisotropic pairing gap by 
the a-model and d-wave pairing. Figure 3 also shows 
^-behavior in these cases obtained by Eilenberger the- 
ory. Even in the gap anisotropy, we find the first order 
transition of orientation, where the ^-dependence is not 
monotonic at low B and Bi st is higher compared to the 
/3-model. Thus, independent of the sources of anisotropy 
(Fermi velocity or gap), there is always the first order 
orientational transition in the similar successive VL tran- 
sition. Starting with 8 — 60° at B c \, the a-direction of 
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Fig. 5. (Color online) Zero energy LDOS N(r, E = 0)/Nq for the 
square VL with a \\ (100) at B/Bgo = 1-2 in the d x 2_ y 2-wave 
case (a) and the /3-model with /3 = 0.5 (b). 



VL coincides with the gap minimum (a-model) or with 
the Fermi velocity minimum (/3-model) . This orientation 
is changed via a first order transition to the orientation 
where the a-direction is rotated by 45°. However, in the 
gap anisotropy case, the reentrant transition from square 
to hexagonal VL at high fields does not occur. We note 
that in these anisotropic gap cases, at low T the non- 
local theory does not work to reproduce the 6*-behaviors 
of Eilenberger calculation. 14 ) 

It is also interesting to notice simple geometry that the 
open angle 8 m [ n — 2tan _1 (|) = 53° when square tiles 
with same size are closely packed as shown in lower inset 
of Fig. 3. In the VL orientation a || (110), as seen from 
Fig. 3, the theoretical minimum of 8 indicates 6> m j n ~ 
53° for both anisotropy cases, which is roughly obeyed 
by the experimental data shown in Fig. 1. In the other 
orientation a \\ (100), 8 -> 90° at B -> B 9Q . This is also 
understandable by the packing of square tiles in different 
way, as shown in upper inset of Fig. 3. 

It is instructive to examine the electronic structures of 
various VL's to understand the origins of the successive 
transition. First, we discuss stable orientation of hexag- 
onal VL at low B in the d a .2_ y 2-wave pairing. In Fig. 4 
we display the typical LDOS at zero energy (i.e., Fermi 
level) for two orientations. In stable orientation a || (110) 
in (a), the zero energy LDOS is well connected between 
nearest neighbor (NN) vortices along a, and between 
next NN vortices. Here a is parallel to node direction 
of d-wave pairing, connected by the A-type bond in in- 
set of Fig. 1. These interconnections effectively lower the 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



kinetic energy of quasi-particles, leading to a stable VL 
symmetry and orientation. In contrast, in the VL of un- 
stable orientation in (b), the interconnections are not well 
organized, clearly demonstrating it less favorable orienta- 
tion and VL symmetry energetically. We note that, even 
in the stable hexagonal VL [Fig. 4(a)], four B type bonds 
defined in inset of Fig. 1 among the six NN vortex bonds 
are not favorable directions for interconnections. Thus 
those are frustrated, which ultimately leads to further 
successive VL transitions in higher fields. 

The concept of interconnections between NN vortices 
via zero energy LDOS continues to be useful for square 
VL at high fields. In d-wave pairing as shown in Fig. 
5(a), interconnections are highly well organized to stabi- 
lize it. In particular not only all NN vortex connections 
are tightly bound, but also second, third neighbor vor- 
tices are also connected. This high field stable VL con- 
figuration can not be continuously transformed from low 
field stable hexagonal VL in Fig. 4(a)just by changing 6. 
Therefore there must exist a first order phase transition 
in an intermediate field region to rotate the orientation 
by 45°. This intuitive explanation is basically correct for 
a superconductor with the four-fold anisotropy coming 
from the Fermi velocity anisotropy. 

In the /3-model shown in Fig. 5(b) in the stable square 
VL, the zero energy LDOS extends broadly to neigh- 
bor vortices. This weak connection between vortices ul- 
timately leads to the instability for this square VL in 
the /3-model towards a hexagonal VL in higher fields as 
shown in Fig. 2. This may be one of the reasons for the 
reentrance phenomenon in the /3-model. Note that in the 
a-model there is no reentrance up to B c2 and the square 
VL is most stable at high fields. 

Let us examine the experimental data in light of the 
present calculations. VsSi 4,5 ) whose data are shown by 
circle symbols in Fig. 1 and Nb 3 Sn 10 ) are cubic crys- 
tals with A-15 structure and known to be an s-wave su- 
perconductor with an isotropic gap. The main four-fold 
anisotropy comes from the Fermi velocity, thus an exam- 
ple of the /3-model. YNi 2 B 2 C 6 ~ 8 ' (the squares in Fig.l) 
has tetragonal crystal and is known to be d xy -\ike gap 
symmetry, 18 ) thus an example of the a- model albeit the 
Fermi velocity anisotropy may be also working. In fact, 
the successive transition precisely follows our calculation, 
ending up with the square VL whose NN is directed to 
(100). 6 ~ 8 ) There is no indication for the reentrance tran- 
sition up to B c2 . CeCoIn 5 9 ) (the triangles in Fig.l) has a 
d x 2_ y 2 gap symmetry. The successive transition, includ- 
ing first order, ends up with the square VL which ori- 
ents along (110) as expected. Towards B c2 it exhibits the 
reentrant transitions to the hexagonal VL with the same 
orientation near B c \. Similar reentrant VL transition is 
also found in TmNi 2 B 2 C. 11 ' This intriguing reentrance 
phenomenon, which is not covered here, belonging to a 
future problem because we need to consider the Pauli 
paramagnetic effect. 19,20 ' 

In conclusion, universal behaviors of the first order 
transition in the VL morphology from hexagonal to 
square lattice have been studied by the non-local Lon- 



don theory and Eilenberger theory. We find that non- 
local London theory is accurate for the superconductor 
whose anisotropy comes from the Fermi velocity, while 
it is uncontrollable for the gap anisotropy case, which 
is covered by Eilenberger theory. We have explained the 
successive VL transition with first order one observed in 
various four-fold symmetric superconductors. Since the 
origin of this phenomenon is due to either gap anisotropy 
or Fermi velocity anisotropy, which are present in any 
superconductors, it is desirable to carefully perform ex- 
periments of neutron diffraction, ^iSR, NMR or STM 
to observe this generic phenomenon. In particular it is 
interesting to see it in high T c cuprates, Sr 2 Ru04 and 
TmNi 2 B 2 C for H/ /(001), which are known to have four- 
fold gap anisotropy, yet so far this phenomenon has not 
been reported although square VL's are observed. 
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